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Summary. We consider an individual-based two-dimensional spatial model with 
nearest-neighbor preemptive competition to study front propagation between an 
invader and a resident species. In particular, we investigate the asymptotic front 
velocity and compare it with mean-field predictions. 



1.1 Introduction and Model 

The dynamics of propagating fronts are fundamental in the study of the 
spread of advantageous alleles, species [1], or opinions [2]. Most notably, 
Fisher [3] and Kolmogorov et al. [4] first addressed the velocity character- 
istics of a simple front by way of a reaction-diffusion equation [1], which 
served as a one-dimensional model for the spread of a favorable gene. Our 
study of front propagation envisions introduction of an advantageous allele 
or a competitively superior species through mutation within [5, 6] or geo- 
graphic dispersal to [7, 8] a resident population, respectively. Introductions 
occur rarely, and stochastically in both space and time. We have shown [5-8] 
that the time evolution of such systems can be well described within the 
framework of homogeneous nucleation and growth. In particular, in two di- 
mensions, for sufficiently large systems, the typical time scale (lifetime) until 
competitive exclusion of the weaker allele or species scales as T~(iu 2 ) -1 / 3 , 
where / is the stochastic nucleation rate per unit area of the successful clus- 
ters of the better competitor and v is the asymptotic radial velocity of the 
corresponding circular fronts. It is, thus, clear that the full understanding of 
the dependence of the lifetime on the local rates of the systems requires that 
of the velocity of the front separating the two alleles or species. Therefore, we 
focus on the velocity characteristics of a two-dimensional two-species model 
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of invasion with preemptive competition, where the invading species has a 
reproductive advantage over the residents. 

Here we study the velocity of the invading fronts of both planar and circu- 
lar shapes on L x x L y two-dimensional lattices as a function of the reproduc- 
tion and mortality rates of each species, and compare them with mean-field 
predictions. Each lattice site can be empty or occupied by the resident or 
the invader. A lattice site represents the minimal amount of locally available 
resources which can sustain an individual. Competition for resources is pre- 
emptive, and therefore an individual site cannot be taken by either species 
until its occupant's mortality makes it available. Preemptive competition is 
typical for plant species competing for common limiting resources [9-13]. The 
local occupation number at site x is rij(x) = 0, 1, i = 1,2, representing the 
number of resident and invader species, respectively. Due to the excluded 
volume constraint, ni(x)ri2(x) = 0, since two species cannot simultaneously 
occupy the same site. A species can colonize open sites through local clonal 
propagation only. An individual of either species occupying site x may repro- 
duce only if one or more of its neighboring sites is empty (here we consider 
nearest-neighbor interactions only). 

Our unit time is one Monte Carlo step per site (MCSS) during which L x L y 
sites are chosen randomly. The occupancy (local configuration) of a chosen 
site is updated based on the following transition rates. When a site is empty, it 
can become occupied by species i of a neighboring site, at rate a^r^x), where 
on is the individual-level reproduction rate and rji(x) = (1/4) X)x'enn(x) n i( x ') 
is the density of species i around site x; nn(x) is the set of nearest neighbors 
of site x. If a site is occupied by an individual, it can die at rate fi (regardless 
of the species). We summarize the local transition rules for an arbitrary site 
x as 

ai ^ x) l, Q2 ^i x) 2, 1-!U0, 2-^0, (1.1) 

where 0,1,2 indicates whether the site is empty, or occupied by the resident, 
or invader species, respectively. We study the regime where [i < ct\ < a 2 , so 
that competition between the two alleles drives the dynamics and the invading 
species has an individual-level reproductive advantage over the resident. 

For planar fronts, we consider an L x xL y lattice with periodic boundary 
conditions along the y direction. The direction of propagation is along the 
x direction by virtue of the initial condition; we set a flat front separating 
the invader species from the residents; for simplicity, to the left (right) of 
the front, all sites are occupied by the invaders (residents). As the simulation 
begins, a number of individuals die in a few time steps, and in both domains, 
away from the front, the densities quickly relax to their "quasi-equilibrium" 
value where the clonal propagation is balanced by mortality. The competition 
between the two species, hence, takes place in the interfacial region. Through- 
out the simulation, we keep track of the location of the invading front, by 
defining the edge as the right-most location of an individual of species 2, for 
each row of our lattice. The average position of the front is then recorded for 
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each time step, from which one can extract the front velocity. We also studied 
the velocity of circular fronts on an LxL lattice, with an initial condition of 
a sufficiently large central cluster of the invading species (with radius slightly 
larger than a critical radius [6,8]), with all other sites occupied by the resi- 
dent species. We then keep track of the time-dependent global density of the 
invaders, from which extracted the average radial velocity of the growing cir- 
cular cluster. Before discussing our simulation results, we first consider those 
obtained from the mean-field equations of motion. 

1.2 Mean-Field Equations and Propagation into an 
Unstable State 

From the above transition rates Eq. (1.1) and the underlying master equation, 
neglecting correlations between the occupation numbers at different sites, for 
the ensemble-averaged local densities p,(x, t)=(rii(x, t)) one obtains 

Pi(x,t + 1) - Pi(x,t) = [1 -pi(x,i) -p 2 (x,t)] ft( x '^) 

x'enn(x) 

-/iA(x,t), (1-2) 

i = 1,2. For further insight we take a naive continuum limit of Eq. (1.2), 
yielding 

d t pi(x,t) = ^ [1 - pi(x,i) - p 2 (x,i)] V 2 pi(x,i) 

+ a t [1 -pi(x,i) - p 2 (x,t)]pi(x,t) - p,pi(x,t) . (1.3) 

Spatially homogeneous solutions of these equations, (p*,p 2 ), are (0,0), (1 — 
p,/ai,0), and (0, 1 — p/a 2 ). For p < ot\ < a 2 , only the (0, 1 — p/a 2 ) fixed 
point is stable. Thus, the motion of the invading front amounts to propaga- 
tion into an unstable state [3,4,14-16], a phenomenon that has generated 
a vast amount of literature [17] since the original papers by Fisher [3] and 
Kolmogorov et al. [4] , with applications ranging from reaction-diffusion sys- 
tems [18,19], population dynamics [1], epidemics [20] or opinion formation 
in social systems [2]. At the level of the mean-field equations, the front is 
"pulled" by its leading edge, and for sufficiently sharp initial profiles, the 
asymptotic velocity v is determined by this infinitesimally small density 
of invaders intruding into the linearly unstable, resident-dominated regime. 
Performing standard analysis on Eqs. (1.3), one obtains the velocity of the 
"marginally" stable invading fronts [1,15-17] 

v* = — Ja 2 (a 2 - ai). (1.4) 
ai 

Thus, for small differences in the local reproduction rates, v*^(a 2 — ai) e 
with 9 = 1/2. Equation (1.4) reproduces the velocity obtained by numeri- 
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Fig. 1.1. (a) Front velocity obtained by numerical iteration of the mean- field equa- 
tions (1.2), for fixed Q2=0.70, as a function of the difference of propagation rates, 
Q2— Qi, for three different values of /i. The dashed curves are the analytic veloci- 
ties of the "marginally stable" fronts, given by Eq. (1.4); the dashed straight line 
segment corresponds to the slope #=0.5. (b) Front velocity from Monte Carlo sim- 
ulations for fixed 02=0.70 as in (a), with L x = 1000, L y = 100, for planar, and 
L x =L y = 1000 for circular fronts, for three values of \x. The dashed straight line 
segment corresponds to the effective power law with an exponent #«0.62 for small 
differences between the local reproduction rates. 

cally iterating the discrete-time discrete-space continuous-density mean-field 
equations (1.2) remarkable well, as shown in Fig. 1.1(a). 

1.3 Monte Carlo Results and Discussion 

We performed dynamic Monte Carlo simulations using the local rates given 
by Eq. (1.1). In the case of planar fronts, we found that the front velocity is 
much smaller than that of the mean- field approximation [Fig. 1.1(b)]. Further, 
for small differences between the local reproduction rates, v*^(ct2 — ai) e 
with 8 w 0.62, an exponent significantly differing from the mean-field scaling 
[Fig. 1.1]. Results from the simulations of the propagation of circular fronts 
closely match those of the planar fronts. A recent study [20] has found a 
similar behavior in a discrete two-dimensional stochastic epidemic model. 

The discreteness of the individuals (or equivalently, effective density cut- 
offs in a continuum description) [21-23] and noise [24,25] have been shown to 
produce velocity characteristic drastically different from those of the mean- 
field equations. More precisely, advancing fronts in stochastic individual- 
based or particle models, which in the mean-field limit converge to a pulled 
front behavior, are instead "pushed" [17]. That is, the front velocity is de- 
termined by the full non-linearity of the frontal region, as opposed to the 
infmitesimally small leading edge [17]. Our model provides an example for 
this generic behavior. 
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